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l^r manUe be more viscous than the upper mantle in order to 
^uce the required positive geoid anomalies. This has already 
been show to be true for the Earth, where the observed geoid highs 
ovCT regions of mantle upwelling and regions of mantle down- 
nf ^ explained by the presence of a strong lower mantle 

U .12], The large positive GTRs and the presence of large shield 
volcanos m certain highland regions on Venus, such as Beta Regio 
Md EisUa Regio, arc best explained as areas of mantle upwelling 
[5.13.14], The regimeof rapid crustal flow predicts crustal thinning 
over the upwelling. However, the extensive partial melt and ensuing 
volcanism expected over such regions of mantle may outweigh the 
effects of crustal thinning on the surface topography and thus also 
yield postive GTRs [15], 
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The first significant dataset to provide local measures of venu- 
sian gravity field variations was that acquired from the Pioneer 
Venus Oibiter (PVO) during the 1 979-1 98 1 period. These observa- 
tions were S-band Doppler radio signals from the orbiting space- 
craft received at Earth-based tracking stations. Early reductions of 
Aese data were performed using two quite different techniques. 
Estimates of theclassical spherical harmonics were made to various 
degrees and orders up to 1 0 [ 1 .2.3]. At that time, solutions of much 
higher degree and order were very difficult due to computer 
limitations. These reducuons. because of low degree and order, 
revealed only the most prominent features with poor spatial resolu- 
tion and very reduced peak amplitudes. 

Another reduction technique was the line-of-sight acceleration 
mapping that had been used successfully for the Moon and Mars 
This approach provided much more detail and revealed the high 
correlauon of gravity with topography [45.6]. However, this tech- 
nique does not produce a global field as do the spherical harmonics. 

It provided a mapping of features from approximately 50°N to 25 °S 
latitude for 360“ of longitude. Other shortcomings were that the 
accelerations were at spacecraft altitude rather than at the surface 
and were not vertical accelerations; however, the reductions were 
quick and cheaply accomplished. Other efforts to analyze these data 
included local area reductions, where surface masses were esti- 
mated [7,8,9]. 

The computer revolution over the past 10 years has allowed new 
reductions with spherical harmonics. New fields up to degree and 
order fifty (2600 parameters) have been made [10 11 12] These 
fields now provide the best representation for any serioui geophysi- 
cist doing quantitativemodeling. There is now vertical gravity at the 
surface from a global model that carries all the requirements of 
dynamical consistency. There is one sizeable concern in that the 


^olution over the entire planet is not uniform. This is due to the 
Pionew orbit, which had a high eccentricity, causing the high 
lauiude regions of Venus to be poorly resolved. 

The Magellan (MGN) spacecraft, which went into orbit about 
V«us m August 1990. has returned Doppler data for gravity field 
^ucUOT. However, because the high gain antenna was pointed at 
Venus for SAR mapping, no gravity data were acquired until the 
antenna wm pointed back to Earth. This occurred at spacecraft 
altitudes higher than 2500 km, gready reducing local gravity 
s^itivity. MGN has an eccentricity much smaller than PVO. so 
mere is new information in the polar regions. Present reductions 
mclude two MGN circulations (486 days), which reduce uncertain- 
lies and produce somewhat better resolution. 

During March. April, and May 1992 new low-altitude data have 
b^n acqu^ ^ both PVO and MGN. PVO periapsis latitude has 
ch^ged 27». from 16“N to 11”S. These data will provide better 
*°'**“" hemisphere, particularly over Artemis. 
Die MGN mission now acquires periapsis gravity data for one orbit 

’ ^“'‘' 80 cs SAR mapping for one orbil/day). Since 
MGN has an X-band radio signal, the data quality is a factor of 10 
tetter dian PVO. Only a small block of MGN data was acquired 
before its penapsis w«ml into occultation May 1 6. Solar conjunction 
and penapsis occultation has also occurred for PVO 

In September of 1992 MGN periapsis will exit occultation and 
US penapsis altitude wiU be lowered to approximately 170 km. 
Penapsis will be visible from Earth for a complete 360° longitude 
average period (243 days). TTiis should be an excellent dataset 
havmg low X-band data noise that in turn can be combined with the 
PVO dataset 

In De^ber 1992 PVO wiU exit periapsis occultation and low- 
^titude data (-150 km) m the southern hemisphere will be acquired 
for about one month before PVO is lost due to the lack of fuel to 
maneuver to safe altitudes. 

In May 1993 there remains the possibility of aerobraking MGN 
mto acircular orbit, thus allowing global uniform resoluuon gravity 
data to be acquired. One hopes that NASA has enough foresight to 
keep Magellan alive so this is a reality. It is anticipated that if this 
IS done, harmonic solutions to degree and order 60-70 (5000 
parameters) will be produced. One could then compare similar 

features globally, resolve coronae and test many interior structure 
models. 
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DIFFERENT TYPES OF SMALL VOLCANOS ON VENUS 
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Ii^Utute, Russian Academy of Science. 1 17975 Moscow, Russia 
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One of the studies of volcanic activity on Venus is the compari- 
son of that with the analogous volcanic activity on Earth. The 
pre iminary report of such a comparison and description of a small 
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Fig. 1. Sketchy cross-acclion of flai-iop volcanos: I -shield edifice; 
2-cinder flat-lop supersimauTC. 
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Fig. Zo Distribution of pit diameters relative to edifice diameters: (a) usual 
small shields; (h) flat top volcanos (pit diameter relative to shield diameter); 
(c) flat-top volcanos (pit diameter relative to flat-top part diameter). 


Fig. 3. Comparison of cralcr/cdifioe diameter ratio histograms for small 
Venusian volcanos (bonom) and possible terrestrial analogs. 


cluster of small venusian volcanos is rcprcsenled in detail in this 
paper. 

The cluster in question of small volcanos on Venus within Sedna 
Planitia (43®N. 334°E) is localized and compact. The total number 
of small volcanos recognized within the cluster is about 160 
edifices; most of them have a small top pit. Some types of volcanic 
edifices within the cluster are recognized. The major type of 
volcanic edifices is small shield volcanos of Iceland type [1 }. Small 
shield average basal diameter is 3750 m. Volcanos with a flat top 
outlinedby a sleep flank (Fig. l)arc also observed. There are ISflat- 
top edifices 9. 4% of the total number) in the cluster. The basal 
diameter of those varies from ISOOmto lO.lOOmand average basal 
diameter is 6450 m, i.c., they are sufficiently larger on the average 
than usual shield edifices. Another difference from regular small 
shield volcanos is the dependence observed between the crater 
diameter and the edifice diameter (Figs . 2a, b). The right dependence 
between crater diameter and flat-top edifice diameter is also ob- 
served (Fig. 2c). The crater diameter to edifice basal diameter ratio 


is similar to that for shield edifices and iciTCStrial small shields (Fig. 
3) [2]. The crater diameter to flat- top part basal diameter is similar 
to that for iciresCrial cinder cones (Fig. 3). If the slope angle of the 
flat-top flank corresponds to the angle of natural slip (of rest) of 
volcanic cinder material on Earth (it is about 30® [3]) the flat-top 
height can be defined, which varies from 1 30 to 260 m for different 
size edifices. 

Thus, at least there are two different morphological types of 
small volcanos within the cluster. The difference may be caused by 
both eruption conditions including volatile content and different 
magma compound. It is suggested that venusian flat-top volcanos 
may have formed in two stages. The first is the formation of usual 
shield edifice due to effusive eruption and the second is the 
formation of the flat cinder superstructure on the shield top (Fig. 1 ) 
due to subsequently more explosive eruption. 

It should be emphasized that venusian small shields are also 
similar to terrestrial analogs in size (1 ]. The minimal basal diameter 
of a venusian small shield obtained by more than 600 measure- 
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menls of volcanos with accuracy down to pixel is within the 
800— 900-m range. The smallest terrestrial shield edifice formed by 
one enq)tion also has a basal diameter of about 800-900 m [3J. 

The different types of volcanos can often be observed within the 
similar volcanic cluster on Earth. There are a variety of volcanic 
structures in similar clusters of volcanos in Saudi Arabia. The 
volcanic cluster and volcanos have been described in detail by C amp 
and Roobol [4]. The scoria cones arc 0.5 to 1 .3 km in diameter and 
rise as much as 100 m above the surrounding surface. They arc 
typically steep and have open craters with a diameter of 1 00— 400 m. 
Scoria cones arc the largest group of Arabian volcanos with numbers 
up to 75% [4] . The shield volcanos arc 5-10 km in diameter and rise 
as much as 1 150 m above the surrounding surface. The crater/cone 
ratio is smaller, which is typical for shield volcanos. This group 
reaches 12% of the total number of volcanos (4J. 
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CONSTRAINTS ON CklJSTAX RHEOLOGY AND AGE OF 
DEFORMATION FROM MODELS OF GRAVITATIONAL 
SPREADING IN ISHTAR TERRA, VENUS. Suzanne E. 
Smrekar and Scan C. Solomon, Department of Earth, Atmospheric, 
and Planetary Sciences, Massachusetts Institute of Technology, 
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Introduction: Gravitational spreading is expected to lead to 

rapid relaxation of high relief due to the high surface temperature 
and associated weak crust on Venus [1,2]. In this study, we use new 
Magellan radar and altimetry data to determine the extent of 
gravitational relaxation in Ishtar Terra, which contains the highest 
relief on Venus as well as areas of extremely high topographic slope. 
Within Ishtar Terra the only mountain belts found on Venus, Akna, 
Danu, Freyja, and Maxwell Montes, nearly encircle the smooth, 
high (3-4 km) plateau of Lakshmi Planum. Finite-element models 
of this process give expected timescales for relaxation of relief and 
failure at the surface. From these modeling results we attempt to 
constrain the strength of the crust and timescales of deformation in 
Ishtar Terra. Below we discuss observational evidence for gravita- 
tional spreading in Ishtar Terra, results from the finite -element 
modeling, independent age constraints, and implications for the 
rheology and timing of deformation. 

Observations: Magellan data have revealed abundant evi- 

dence for localized extension throughout Ishtar Terra [3,4]. Many of 
the observed exlensional features are oriented perpendicular to the 
downslope direction, implying that they formed as a result of 
gravitational spreading. In some of the mountainous areas, exten- 
sional faults occur parallel to the apparent direction of shortening. 
By analogy with the Himalayas on Earth [5], we interpret such 
extension to indicate graviulional spreading during convergence. 
Areas of possible gravitational spreading are observed in each of the 
four mountain belts and along both the southern and northeastern 
margins of Lakshmi Planum. Horizontal strain in these regions is 
estimated to be 1-20%; this value may be an underestimate if blocks 
are rotated along the normal faults. The observational evidence for 
gravitational spreading in IshtarTerra is discussed more fully in [4]. 

Models: We use the finite-element algorithm TECTON [6] to 

model the evolution of gravitational spreading in a vertical section 
of the crust near the margin of a plateau or the edge of a broad 
mountain belt Each model includes an elevated plateau 50 km in 


width, a sloped margin, and a lowlands region 80 km in width. The 
code employs a depth-dependent, viscoelastic rheology with non- 
linear stress dependence and exponential temperature dependence. 
We adopt a Young’s modulus of 6 x 10^® Pa, a Poisson’s ratio of 
0.25, and adiabase flow law [7]. Although a diabase composition is 
most similar to the composition determined at Soviet lander sites 
[8], the flow law determined for diabase may underestimate the 
strength because of partial melting of samples during tfic experi- 
ment [7,9]. For this reason, we also employ a webstcritc flow law 
( 10] as an approximate i^rper bound on crustal strength. The surface 
temperature is 740 K and increases linearly with depth. Each row of 
elements in the grid has the same viscosity, which is equivalent to 
assuming that the temp>erature is constant along both the bottom and 
lop of the grid. This results in a somewhat higher thermal gradient 
in the plains than in the plateau. The boundary conditions are zero 
vertical and horizontal velocity on the bottom of the grid, a free 
surface top boundary, and zero horizontal velocity and flee vertical 
slip al the sides. The bottom boundary condition approximates an 
upper mantle layer that is much stronger than the lower crust. 
Details of the models are given in [4]. Brittle failure is evaluated 
using a Mohr-Coulomb criterion. The timing of predicted brittle 
failure and relaxation of rcUef are found for ranges of plateau height 
(1—5 km), plateau margin slope (1®— 30®), crustal thickness (10-30 
km), and thermal gradient (5-25 K/km). Slopes of l®-30® and 
plateau or mountain bell heights of 1—6 km arc observed in Ishtar 
Terra, Crustal thickness on Venus is predicted to be 10-30 km 
[11,12]; the average thermal gradient is expected to be KJ-25 K/km 
[13,14]. 

The rate of deformation is largely controlled by the effective 
viscosity, and thus the temperature and stress, al the base of the 
crust. The greatest deformation occurs at the base of the crust, where 
the viscosity is lowest. Temperature al the base of the crust increases 
with crustal thickness (for a fixed thermal gradient), and deviatoric 
stress near the base of the crust increases with both crustal thickness 
and plateau elevation. The nonlinear change in predicted failure and 
relaxation lime (see below) is due both to the exponential depen- 
dence on temperature and to the pK>wer law dependence on deviatoric 
stress of the effective viscosity. 

Failure is predicted to occur when the stresses exceed the Mohr- 
Coulomb criterion. The earliest failure is located on the plateau. In 
most cases, failure on the slope follows shortly after failure on the 
plateau; thrust faulting in the plains is predicted typically only for 
a crustal thickness of 30 km. Figure 1 shows the lime from the start 
of spreading until the first failure on the plateau as a function of 
plateau height and crustal thickness, for a margin slope of 3® and a 
thermal gradient of 15 K/km. Models with both webslerilc and 
diabase flow laws are plotted. No times are shown for models that 
do not fail within 1 Ga. Cases with a diabase rheology, a large 
plateau height, and a large crustal thickness fail almost immedi- 
ately, A plateau height of 1 km and a crustal thicknesses of 10 km 
can lake up to tens of millions of years to fail. With a websterite flow 
law, no failure occurs unless the crust is a least 20 km thick and the 
plateau elevation is 1— 3 km. Changing the thermal gradient by 10 K/ 
km has a large effect on failure limes, with a higher thermal gradient 
decreasing the failure time by a factor of lO^lO^, and a lower 
thermal gradient increasing times by a similar factor. When failure 
is first predicted, the horizontal surface strain is approximately 
0.05%. This strain is unlikely to be recognized as normal faulting in 
Magellan radar images; a strain of 1%, which accumulates once 
significant relaxation of the topography begins, is probably a more 
reasonable value to compare with observations. As faulting is not 
explicitly modeled, there is some uncertainty in the interpretation of 
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